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In this work, we have studied the Brans-Dicke (BD) cosmology in anisotropic models. We present 
three dimensional dynamical system describing the evolution of anisotropic models containing per- 
1 feet fluid and BD scalar field with self-interacting potential. The relevant equations have been 

transformed into the dynamical system. The critical points and the corresponding eigen values have 
^ ■ been found in radiation, dust, dark energy, ACDM and phantom phases of the universe. The natures 

' and the stability around the critical points have also been investigated. 
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I. INTRODUCTION 

A recent renewal of interest in Brans-Dicke (BD) theory [1] can be traced to the discovery by La and 
Steinhardt [2] that the use of BD theory in place of general relativity can ameliorate the exit problem of 
£^ n . inflationary cosmology. This is possible because the interaction of the BD scalar field with the metric slows the 
expansion from exponential to power-law. The BD theory contains only one dimensionless parameter, u> and 
the effective gravitational constant is inversely proportional to the scalar field </>. BD theory has been proved to 
be very effective regarding the recent study of cosmic acceleration [3] . This theory yields the correct Newtonian 
weak-field limit, but solar system measurements of post-Newtonian corrections require u) > 500 [4]. In the limit 
uj — > oo, the field 4> becomes a constant [5] and we recover Einstein gravity. This theory has very effectively 
O ' solved the problems of inflation and the early and the late time behaviour of the Universe. N. Banerjee and 
1—1 1 D. Pavon [3] have shown that BD scalar tensor theory can potentially solve the quintessence problem. The 
generalized BD theory [6] is an extension of the original BD theory with a time dependent coupling function 
bj. In Generalized BD theory, the BD parameter a; is a function of the scalar field <f>. This has led to more 
general scalar-tensor gravity [7] being considered with a self-interacting potential [8, 9]. Modified BD theory 
£f~) , with a self-interacting potential have also been introduced in this regard. Bertolami and Martins [10] have used 
<3\ • this theory to present an accelerated Universe for spatially flat model. They have obtained the solution for 
accelerated expansion with a potential (j> 2 and large |cj|, although they have not considered the positive energy 
£SJ ■ conditions for the matter and scalar field. 

t-h ; 

Dynamical system theory [11] has been applied with great success in cosmology and astrophysics within the 
context of general relativity. This theory are used to describe the behaviour of complex dynamical systems 
, usually by constructing differential equations. This theory deals with a long term qualitative behaviour of 
the formed differential equations. It does not concentrate to find the precise solutions of the system but 
r> | provide answers like- whether the system is stable for long time and whether the stability depend on the 
initial conditions. Besides the other scientific fields this theory is now become widely useful in the research of 
cosmology. Its range of applicability has been enlarged by considering alternate theories of gravity, such as 
those which do not obey the principle of minimal coupling [12], Kaluza-Klein [13] and scalar-tensor theories 
[14]. There are several works on qualitative analysis in dynamical system of FRW cosmology in BD gravity [9, 
15]. There are some literature where the authors have studied dynamical evolution of some models where the 
universe is filled with dark energy on an interaction has occurred between dark matter or barotropic matter 
and some dark energies like Generalized Chaplygin Gas [16], New Generalized Chaplygin Gas[17], ELKO non 
standard spinor dark energy [18], DBI essence [19], non- minimally coupled scalar field [20], k-essence [21] etc. 
There are some other approach of dynamical system analysis to diffrent model the universe on the frame work 
of loop quantum gravity and non linear electro dynamics [22, 23]. Some authors [9, 24] have shown that the 
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field equations for FRW cosmological models could be reduced to a two dimensional dynamical system for any 
reasonable perfect fluid matter source in BD theory and in the presence of scalar potential, the cosmological 
dynamical system turns out to be generally three dimensional, apart from radiation dominated universes in 
which the system once more reduces to two dimensions. 

Here we have considered an anisotropic model [25, 26] containing perfect fluid and BD scalar field with 
self-interacting potential. The starting point of our work is to reduce the field equations to a three dimensional 
autonomous dynamical system. From the knowledge of mathematical features of the system, we have found 
the critical points and the exact solutions for the system of equations. Some of the solutions are valid for fluids 
satisfying the equation of state and for particular values of 7. The stability around the critical points have 
been investigated. 



II. BASIC EQUATIONS 

The self- interacting Brans-Dicke (BD) theory is described by the action: (choosing 8nGo = c = 1) 



d X\/—g 



<f>R--<f>> a 4>, a -V(<t>)+C r . 



(1) 



where V{<j>) is the self-interacting potential for the BD scalar field <f> and ui is the BD parameter. The matter 
content of the Universe is composed of perfect fluid, 



Tpv = {p + p)u^u v + p g^ 
where u v = — 1 and p, p are respectively energy density and isotropic pressure. 
From the Lagrangian density (1) we obtain the field equations 



+ \ - <v □ <j>] - ^£r<w + - T < 



2(f> 



and 



• dV{<t>) 



- 21/(0) 



3 + 2oj 3 + 2oj 
where T = T^g^. 

We consider homogeneous and anisotropic space-time model described by the line element 



(2) 



(3) 



(4) 



ds z = -dt z + a z dx z + b z dtt z k 
where a and b arc functions of time t alone : we note that 



(5) 



dn 2 k 



dy 2 + dz 2 , when k = (Bianchi I model) 

when k = +1 (Kantowaski-Sachs model) 



+ sinhr9d<j) , when k = — 1 (Bianchi III model) 



Here k is the curvature index of the corresponding 2-space, so that the above three types are described by 
Thorne [19] as flat, closed and open respectively. 

Now, in BD theory, the Einstein's field equations for the above space-time symmetry are 



a b 

- +2 l 

a b 



1 



(3 + 2u)(t> 



[(2 + u)p + 3(1 + u)p] - - ^ 



4> 2 j> v{<j>) 



2<j) 



(6) 
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^_ , 2 - - = - - — - ( - +2-1 ^ 4- 4- 

b 2 ab 4> b 2 \a b) 4> 2(f) 2 2(j) 



(7) 



and the wave equation for the BD scalar field (j) is given by 



i 



3 + 2w 



(8) 



The energy conservation equation is 



(9) 



Here we consider the Universe to be filled with barotropic fluid with EOS 



p=( 7 -l)p, 0< 7 <2 
The conservation equation (9) yields the solution for p as 



(10) 



p = po (ab 2 ) 



2\-7 



where, po is an integration constant. 



(11) 



III. DYNAMICAL SYSTEM ANALYSIS 

In this section, we shall define some variables such that the above equations can be transformed to first order 
differential equations and next we investigate the dynamical system analysis. Let us consider, V((f>) = V^4> n+1 
and define the following variables 



X 



(12) 



1 (a J\ 
Y =-[— + 2- + — 

3 V a b 2<P 



(13) 



where, ' = -j- = (ab 2 ) 1 ^ ^ , -q is the conformal time. If X > then Y must be positive. If X < then Y 
may or may not be positive. 

Using the transformations (12) and (13), the equations (6) and (7) become the following two transformed 
equations (for k — 0): 



X' + 2XY = (ab 



2\2/3 



(4 - 37) p n-l 



(3 + 2u>) (f> 3 + 2lo 



V < 



(14) 



Y' + 2Y 2 = (ab 2 ) 



{5 + 3cu - (2 + oj)j} p (10 + 6fa;-n) ,„ 
3 + 2cj 4> 6(3 + 2w) 0(p 



(15) 



Also for simplicity, let us choose n = z 2 _ 2 an d defining another transformation 
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Z = V Q <j>(al> 



2N7-2/3 



(16) 



(which is always positive) the above equations reduce to the system of non-linear first order differential 
equations 



X' + 2XY = 



4-37 
3 + 2w 



Vn ( Z \ — 



Vo 

Z 37-2 \V 



(17) 



, 2y2 = {5 + 3u-(2 + co) 1 } V {(3 7 -2)(5 + 3^)-l}l/ f Z_\ ~ 
3 + 2w Z 3(3 + 2w)(37-2) \V , 



(18) 



z 1 = (37 - 2)zr - — — -zi 



(19) 



Now in special cases, we consider radiation, dust, dark energy, ACDM and phantom models and find the 
critical points, eigen values and its natures separately. 

Case (a): Dark Energy: For 7 = 1/3, the system of equations (20) becomes 



X' = -2XY + W ° , + J V(f , 
Z(3 + 2uj) Z 2 (3 + 2w) 

r = _ 2Y 2 + (13 + MVo , (2 + 3c)y 3 



3Z(3 + 2o>) 3Z 2 (3 + 2o;) 



Z' = ^XZ - YZ (20) 

There exists two critical point namely (X t DE , YP E , Zf E ) | i=lj2 = (zpffi, T -^L, 3(5 ~^" ) 1 V '° ) where /3 = 
y ( i6.-i)( 7+ 5^ _ and 1 5 

A /(5-2w){15-4(w-l)w} 16 2 

Characteristic Equation around the critical point (XP E ,Y i DE , ZP E ) (,=1,2 are given by the equation 

£ 3 + A H £ 2 + A 24 £ + A 3t = 0, i = 1,2 (21) 

where 

v/2(16w- l)(7 + 5w) 



6(5 - 2w)i(3 + 2w)V r V±5 - 4(w - l)w 



. (16a;- 1)(32cj 2 - 442o>-449) . 

A 2 i = =F 7 — z — 1,2 

324(5 - 2cj) 4 (3 + 2w) 2 y 2 

A . =T (16^ -1)1(7 + 5c) 

3i 81(5 - 2w)i(3 + 2w) 2 y V2{15 - 4(o> - l)o>} 

The condition for the existence of stable attractor solution for the dynamical system (20) is An > 0, A 2 i > 
and AuA 2 i - A 3l > , for i = 1, 2. 

In our model, the dynamical system (20) is not stable around (XP E , YP E , Zf 5 ') but stable attractor in the 
late time universe around the critical point (XP E , YP E , ZP E ). Corresponding condition becomes as follows: 

(16a; - 1)5(7 + 5w)2 (499 + 442o> - 32o> 2 - 12(5 - 2o>) 3 (3 + 2o>)(16o> - l)V H 

n > (23) 

972\/2(5-2cj) j #(3 + 2u;) 3, 1/ 2 {15 - 4 (a; - l)u>}3 
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In particular, for our choice 7=5, w = 0.5, Vb = 0.5, critical point is (0.679563,1.01934,0.428571) and 
the corresponding eigenvalues are (—4.07738,-1.85528,-0.183411), which shows a stable attractor. Another 
critical point is (—0.679563,-1.01934,0.428571), so the system is unstable around that fixed point as the 
eigenvalues are (4.07738,1.85528,0.183411). The phase space diagram of parameters X(rj), Y(rj), Z(rj) and 
their progressions have been drawn in figures 1 and 2 respectively. 
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FIG. 1: The phase space diagram of parameters X(r]), Y{rf), Z(rj) for 7 = i,w = 0.5, Vo = 0.5. The initial conditions 
chosen are X(0) = 0.1, Y(0) = 0.9, Z(0) = 0.3V/ (green); X(0) = 0.2, F(0) = 1, Z(0) = 0.4V/ (blue); X{0) = 0.3, 
F(0) = 1.1, Z(Q) = 0.5V/ (red); X(0) = 0.4, F(0) = 1.2, Z(0) = 0.6^ (brown). 

FIG. 2: The progression of X(rj), Y(r)), Z(t]) for 7 = \,w = 0.2, Vo = 2 and initial condition is X(0) = 0.4, Y(0) = 0.6, 
Z(0) = F § . 

Case (b): Phantom: For 7 = —1/2, the dynamical equation (20) becomes: 



X' 



Y' 



-2XY 



llVn 



nv 



z 



-2Y 



2(3 + 2cj)Z 7(3 + 2lj)Z\V 
., (12 + 7w)V (37 + 21w)Vb 



2(3 + 2cj)Z 21(3 + 2w)Z 
Z' 



7 11 



(24) 



For a; = — i, there exists critical point (X p ,Y P , Z p ) (superscript "P" stands for phantom) for positive 
potential function Vq. 

In particular, 7 = — i w = —0.5, Vq — 4, critical point is (2.8606,2.24761,1.497) and the cor- 
responding eigenvalues are (—8.99046,-4.46316,-0.0320722), which shows a stable attractor. Another 
critical point is (—2.8606,-2.24761,1.497) system is unstable around that fixed point as the eigenvalues are 
(8.99046,4.46316,0.0320722). The phase space diagram of parameters X(rj),Y(rf), Z(rj) and their progressions 
have been drawn in figures 3 and 4 respectively. 
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FIG. 3: The phase space diagram of parameters X(rj),Y{rf), Z{rf) for 7 = — |,w = —0.5, Vo = 4. The initial conditions 
chosen are X(0) = 2.2, Y(0) = 2.4, Z(0) = l.2v} (green); X(0) = 2.3, Y(0) = 2.6, Z(0) = 1.3V * (blue); X(0) = 2.4, 
y(0) = 2.8, Z(0) = 1.4v/ (red); X(0) = 2.5, Y(0) = 3, Z(Q) = l.5V$ (brown). 

FIG. 4: The progression of X{rj), Y(rj), Z{n) for 7 = -|,w = -0.5, V = 2 and initial condition is X(0) = 0.8, 
F(0) = 0.6, Z(Q) = . 

Case (c): Radiation: For 7 = 4/3, the system of equations (20) becomes 



X' = -2XY 

Y> = -2Y 2 4- £ 4- ( 7 + 5 "W> 
2 3Z(3 + 2w) 

Z' = 2yz (25) 

From the above dynamical system of equations, we can seen that there is only one feasible critical points 
namely (X R , 0, \Z2aVo) ("i?" stands for Radiation) where a — — 3 ( 7 3 | f 5 2 L ^) and the feasible range for w becomes 

— I < w < — |, A^ 7^ is arbitrary real constant. In the feasible choice 7 = |, w = —1.45, Vb = 0.5, only 
critical point becomes (0.1,0,0.912871), eigen values are (1.3512,-1.3512,0) depicts unstable solution. The 
phase space diagram of parameters X(rj), Y(rj), Z(r\) and their progressions have been drawn in figures 5 and 6 
respectively. 

Case (d): ACDM: For 7 = 0, the system of equations (20) becomes 



X' = -2XY + W ° T+ W « , 
Z(3 + 2lj) Z(3 + 2uj) 

Y' = -2Y 



, (5 + 3uj)V (11 + 6u)V£ 



Z(3 + 2w) 6Z(3 + 2w) 
Z' = 2XZ - 2YZ (26) 

From the above dynamical system of equations, we can seen that there is only one feasible critical 
point namely (X A ,Y A ,Z A ) ("A" stands for ACDM) where X A — Y A = - V^+^+ a^ Z A > 

(1— &u)y/ Z A (Z+2bj)) 

is arbitrary real constant and the feasible range for uj becomes — | < u < i. In the feasible choice 
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FIG. 5: The phase space diagram of parameters X(rj), Y(rf), Z{rf) for 7 = |,w = —1.45, Vo = .5. The initial conditions 
chosen are Z(0) = 0.2, Y(0) = 0.4, Z(0) = V * (green); X(0) = 0.4, Y(0) = 0.6, Z(0) = (blue); X(0) = 0.6, 
Y(0) = 0.8, Z(0) = V * (red); X(0) = 0.8, Y(0) = 1, Z(0) = (brown), X(0) = -0.2, Y(0) = 0.4, Z(0) = (green); 
X(0) = -0.4, Y(0) = 0.6, Z(0) = V * (blue); X(0) = -0.6, Y(0) = 0.8, Z(0) = ^ (red); X(0) = -0.8, Y(0) = 1, 
Z(0) = (brown). 

FIG. 6: The progression of X(ri), Y(rj), Z(rf) for 7 = |,w = —1.45, Vo = 1.5 and initial condition is X(0) = 0.4, 
Y(0) = 0.6, Z(0) = v} . 

7 = 0, u) = -0.1, Vb = 0.2, only critical points become (±0.463304, ±0.463304, 20.2) and the corresponding 
eigenvalues are (=Fl. 85322, =p0. 931517, ±0.00490953) and so the critical points are unstable. The phase space 
diagram of parameters X(j]), Y(j]), Z(rj) and their progressions have been drawn in figures 7 and 8 respectively. 

Case (e): Dust: For 7=1, the system of equations (20) becomes 



X' = -2XY 
Y' - 



Vo 



-2Y^ 



Z(3 + 2cj) Z(3 + 2co)V 
,2,^0 Z 2 (4 + 3^) 
Z 3V (3±2u;) 

Z> = ™ + YZ 



(27) 



For the dust dominated universe, there are two critical points: (X- D ,Y i D ,Z[ > ) |i=i,2, where U D" stands for 
dust dominated universe, and given by 



X? 2 = T 



1 i.2 r 



2(13 + 9w)V 3 



(7±4w)*x/3 + 2w (5 + 6uj)i 



(13 + 9w)V 3 



(7 + 4w)s v , 2(3±2w) (5±6w)s 
(21 + 12w)3V ^ 



7 D - 



(5 + 6w)3 



(28) 
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FIG.7 



FIG.8 



FIG. 7: The phase space diagram of parameters X(ri), Y(rj), Z(rj) for 7 = 0,cj = —0.1, VI) = .2. The initial conditions 
chosen are X(0) = 3.2, Y(0) = 3.2, Z(0) = 60.5V/ (green); X(0) = 3.4, Y(0) = 3.4, Z(0) = 60.6V ^ (blue); X(0) = 3.6, 
Y(0) = 3.6, Z(Q) = 60.7V/ (red); X(0) = 3.8, Y(0) = 3.8, Z(0) = 60.8V/ (brown). 

FIG. 8: The progression of X (7?), Yfa), ^fa) for 7 = Q,w = -0.2, V = 2 and initial condition is X(0) = 0.4, F(0) = 0.2, 
Z(0) = V § . 



Feasible region for existence of these critical point is u; > — |. Characteristic roots near the critical point 
(Xi , Y®, Z® ) are two positive and one negative real roots, and around the critical point (X^ , Yf , Z®) are 
two negative and one positive real roots. In the feasible range 7 = 1, ui = —0.6, Vq = 0.2, critical points are 
(±1.28263, ±0.641315, 0.733281), and the corresponding eigen values are (=f3. 17407, ±0.696701, T0.0878872) 
and so they are unstable. The phase space diagram of parameters X(j)), Y(j)), Z(rj) and their progressions have 
been drawn in figures 9 and 10 respectively. 



IV. DISCUSSIONS 



In this work, we have studied the Brans-Dickc (BD) cosmology in anisotropic models. We present three dimen- 
sional dynamical system describing the evolution of anisotropic flat (k = 0) model containing perfect fluid with 
barotropic EoS and BD scalar field with self-interacting potential. Choosing the power law form of the potential 
function V in terms of 4>, and defining the three variables X, Y and Z, the field equations can be transformed 
into the dynamical system. The critical points and the corresponding eigen values have been found in radiation, 
dust, dark energy, ACDM and phantom phases of the universe. The natures and the stability around the 
critical points have also been investigated. For dark energy case, we have considered 7=5, w = 0.5, Vq — 0.5, 
and one critical point is (0.679563,1.01934,0.428571), which is a stable attractor. Another critical point is 
(—0.679563,-1.01934,0.428571), which is unstable around that fixed point. The 3D phase space diagram of 
parameters X, Y, Z have been drawn in figure 1 corresponding to the stable critical point. From figure 2, we see 
that X and Y initially increase and then decrease and Z increases in late stage of the universe. For phantom 
case, we have taken 7 = — |, ui = —0.5, Vo = 4, and one critical point is (2.8606,2.24761,1.497), which is 
a stable attractor. Another critical point is (—2.8606,-2.24761,1.497), which is unstable around that fixed 
point. The 3D phase space diagram of parameters X, Y, Z have been drawn in figure 3 corresponding to the 
stable critical point. From figure 4, we see that X and Y initially increase and then decrease and Z increases in 
late stage of the universe. For radiation case, we have considered 7 = |, u = —1.45, Vo = 0.5, and the only 
critical point becomes (0.1,0,0.912871) which is unstable. The 3D phase space diagram of parameters X, Y, Z 



FIG.9 
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FIG. 9: The phase space diagram of parameters X(rf), Y(rj), Z(rj) for 7 = l,w = —0.6, Vq = .2. The initial conditions 
chosen are X(0) = -1.2, F(0) = -0.2, Z(0) = 1.2V/ (green); X(0) = -1.1, Y(0) = -0.3, Z(0) = 1.3V/ (blue); 
X(0) = -1.0, Y(0) = -0.4, Z(0) = 1.4V/ (red); X(0) = -0.9, Y(0) = -0.5, Z(0) = 1.5V/ (brown). 

FIG. 10: The progression of X(tj), F(t?), #(77) for 7 = 1, w = -0.7, Vo = 2 and initial condition is X(0) = 0.4, Y(0) = 0.6, 
Z(0) = v/ . 

have been drawn in figure 5 corresponding to the critical point. From figure 6, we see that Y and Z increase 
and X decreases. For ACDM case, we have taken, 7 = 0, ui = —0.1, Vo = 0.2, and the critical points become 
(±0.463304, ±0.463304, 20.2) and these are unstable. The 3D phase space diagram of parameters X, Y, Z have 
been drawn in figure 7 corresponding to the critical point. From figure 8, we see that X and Y initially 
increase and then decrease and Z increases in late stage of the universe. For dust case, we have assumed 
7=1, uj = -0.6, V = 0.2, and the possible critical points are (±1.28263, ±0.641315, 0.733281) which are 
unstable. The 3D phase space diagram of parameters X, Y, Z have been drawn in figure 9 corresponding to the 
critical point. From figure 10, we see that Y and Z increase and X decreases. So the anisotropic model of the 
universe in Brans-Dicke theory can be stable for some cases of the fluid distribution in late stage of the evolution. 
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